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We are interested in studying how different force fields interact with gravity, at a 

"fundamental" level, but not at the level of Quantum Field Theory. This is because there 

'—i ■ exists no theory of quantum gravity, and no understanding of "Planck-scale" physics; that 
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where G is Newton's gravitational constant, h is Planck's constant, and c denotes the 
$Li | speed of light. 
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Before discussing our results, we think that it is worthwhile to see what is the difficulty 
in making a theory of quantum gravity. In order to understand why a solution to the 
problem of reconciling gravity and quantum mechanics has been so elusive, we must 
consider the implications of the Heisenberg Uncertainty Principle (HUP) at small distance 
scales. 

Recall that the HUP states that "the more precisely a spatial measurement is made, 
the less precisely the momentum (or the energy) of the system being measured, is known" . 
When the spatial measurement Ax ~ 1CU 13 cm, there are large uncertainties in the energy 
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E. These are realized as "fluctuations" at small distances. Since E = mc 2 , the energy 
of these fluctuations can give rise to "virtual" particles and anti-particles, which arise 
out of the vacuum for short times before annihilating each other, and give rise to a "sea 
of virtual particles". When Ax ps 10~ 33 cm (Planck length), the energy fluctuations are 
quite large over small distances, and general relativistic effects become important. In fact, 
from a theorem of Schoen and Yau [12], singularities must form, and these are believed to 
be black holes. Thus space-time becomes very curved at small distance scales; physicists 
say that space-time is "foamy" (not smooth). This invalidates the usual computational 
techniques of Quantum Field Theory, where small curvature is needed - the calculations 
break down at high energies. This suggests that at high energies, General Relativity, or 
Quantum Mechanics, (or both) must be modified. The unsolved problem is, how is this 
to be done? 

In order to set the background for our discussion, we shall briefly review what we con- 
sider to be some of the most important results concerning the coupling of gravity to other 
fields. The first such result is due to Reissner, and Nordstrom (1918, 1919) whereby they 
coupled gravity to electromagnetism, and this led to the celebrated Reissner-Nordstrom 
solution, about which we shall have more to say in the next sections. In the 1920 's, 
Einstein, de Sitter, Friedmann, and others coupled gravity to perfect fluids, in order to 
study problems in Cosmology, the large scale structure of the Universe. In the 1930's Op- 
penheimer, Tolman, Snyder, Volkoff, Landau, and others also coupled gravity to perfect 
fluids in order to study problems connected with stellar formation and collapse. (All of 
the above results are classical and can be found in most standard textbooks on General 
Relativity; see e.g. [1].) 

In the more modern era, the first important ideas stem from the 1988 paper of R. 
Bartnik and J. McKinnon, [2]. These authors coupled Einstein's equations to an SU(2) 
Yang-Mills field (that is, gravity to the weak nuclear force), and they found, numeri- 
cally, particle-like solutions. In the papers [15, 16], Smoller, Wasserman, and Yau rig- 
orously proved the existence of both particle-like and black-hole solutions for the SU(2) 
Einstein/ Yang-Mills equations. Also in recent years, Smoller and Temple [13, 14] coupled 
gravity to perfect fluids and they constructed a theory of shock-waves in General Rel- 



ativity, with applications to Cosmology and stellar structure. In the present paper, we 
couple gravity to both spinors and electromagnetism; i.e., we couple Einstein's equations 
to both the Dirac equation and Maxwell's equations and we investigate the existence 
of black-hole solutions; (in the works [4, 5], we consider particle-like solutions of these 
equations). We shall show here that in all cases considered, black-hole solutions do not 
exist; for complete details, see [6, 7], and also [8, 9] for related work. 



2 Background. 

In this section we shall give a short discussion of Einstein's equations of General 
Relativity, Maxwell's equations of electromagnetism, and Dirac's equation of relativistic 
Quantum Mechanics. For more details, the reader should consult the standard textbooks, 
e.g. [1, 11]. 

We begin with General Relativity (GR). The subject of GR is based on Einstein's 
three important hypotheses: 

(Ei) The gravitational field is described by the metric gij in 4-d space-time; g^- is assumed 
to be symmetric: g^ = g^, i,j = 0, 1, 2, 3. 

(E 2 ) At each point, the 4x4 symmetric matrix [g^] can be diagonalized as g^ = diag 
(1,-1,-1,-1). 

(E 3 ) The equations of GR should be independent of the choice of coordinate system. 

The hypothesis (Ei) is Einstein's great insight, whereby he "geometrizes" the gravi- 
tational field. (E 2 ) means that there should exist "inertial frames" at each point, so 
that Special Relativity is included in GR, and (E 3 ) implies that the gravitational field 
equations should be tensor equations. 



The metric g^ = gij(x), i,j = 0, . . . , 3, x = (x ,^ 1 ,^ 2 ,^ 3 ), x° = ct, is the metric 
tensor defined on 4-d space-time. Einstein's equations are ten (tensor) equations for the 
metric (gravitational field) and take the form 

Rij - -Rgij = crTij . (2.1) 

The left-hand side Gij = Ry — \ R g^ is the Einstein tensor, and is a geometric object, 
while Tij the energy-momentum tensor, represents the source of the gravitational field 
and encodes the distribution of matter. The word "matter" in general relativity refers to 
everything which can produce a gravitational field. The tensor T^ is required to satisfy 
the relation Tl^ = 0, (it's covariant divergence vanishes), and this in turn expresses 
the laws of conservation of energy and momentum. The quantities which comprise the 
Einstein tensor G^ are given as follows: First from the metric tensor g^, we form the 
Levi-Civita connection F k ,: 



,k = 1 kt ^9ii dg ie d 9ij 
lj 2 g { dx i dx i dx e 



where [g ke ] = [g^e] 1 , and summation convention is employed; namely an index which 

^k 



appears as both a subscript and superscript is to be summed form to 3. Having T^ 



we then construct the Riemann curvature tensor R l qk f 
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Finally, we can explain the terms R^ and R in E^; namely R^ = i?f • is the Ricci tensor, 
and the scalar R = g^Rij, is called the scalar curvature. 

The quantity a is a universal constant defined by 

8ttG 
a - 
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where G is the gravitational constant, and c is the speed of light. 

From these definitions, one immediately sees the enormous complexity of the equa- 
tions (2.1) for the unknown quantities g^. For this reason, one seeks solutions which 



are highly symmetric, and in what follows, we shall only consider solutions which are 
spherically symmetric; i.e. the metric has the form 

ds 2 = T-\r) dt 2 - A-\r) dr 2 - r 2 dti 2 , (2.2) 

where dfl 2 = dfi 2 + sin 2 dtp 2 is the standard metric on the 2-sphere, (r, -Q, tp) are the usual 
polar coordinates, and t denotes time. Such metrics turn out to be quite interesting 
mathematically, in addition to having physical interest. 

Now consider the problem of finding the gravitational field exterior to a ball in IR 3 ; 
that is, there is no matter exterior to the ball. This actually models well our solar system, 
where we think of the Sun as the 3-ball, and we ignore the masses of the planets. In this 
case the Einstein equations become Rij — \ Rgij = 0, and we seek spherically symmetric 
solutions. This problem was solved by K. Schwarzschild in 1916, and its solution is called 
the Schwarzschild metric: 

ds 2 =(l-*HL) Jdt 2 -(l- —Y'dr 2 - dn 2 . 



Here m = GM/c 2 and M is the mass of the 3-ball. Since 2m has dimension of length, it 
is called the Schwarzschild radius. One sees immediately that the Schwarzschild metric 
is singular at r = 2m; namely goo = and gu = oo. For r < 2m, the signs of the metric 
components change: g 00 > and gu < 0. So we must reconsider the physical meaning 
of t and r inside the Schwarzschild radius. 

To investigate the mathematical and physical nature of the Schwarzschild radius, we 
introduce Kruskal coordinates (r,t) — ► (u,v), and we seek a transformed metric of the 
form 

ds 2 = f(u, v) (dv 2 - du 2 ) + r 2 dti 2 . 

This yields the following in the region r > 2m, (see Figure HD 
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Figure 1: Kruskal Coordinates 

with similar expressions in the region r < 2m. 

We can use Kruskal coordinates to study light rays traveling radially inward towards 
the Schwarzschild radius, starting at a point P outside of the Schwarzschild radius; i.e., 
in the region r > 2m; see Figure 1. In u/v coordinates, the light ray leaves at the point 
P in r > 2m, t finite, and travels inward towards the Schwarzschild radius r = 2m, as 
t — ► oo. It crosses the line t = oo into the interior of the Schwarzschild sphere. Incoming 
light is thus in effect, totally absorbed by the Schwarzschild sphere. 

We can also study light emitted from inside the Schwarzschild sphere starting at a 
point Q in the region r < 2m; again see Figure 1. The trajectory starts at some r < 2m 
and finite t, travels through increasing r but decreasing t, and crosses the Schwarzschild 
radius t — — oo to the exterior of the Schwarzschild sphere, where its evolution is normal. 
Thus, light emerging from the Schwarzschild sphere must have been traveling since before 
t = — oo; in effect since before the beginning of time. It is questionable whether such 
light is physically observable. If not, the Schwarzschild sphere has the physical properties 
of a black-hole: it absorbs all light and emits none. 

In general, for a metric of the form (2.2), we define a black-hole solution to be a 
solution of Einstein's equations which satisfies 



A{p) = 0, 



and 



A(r) > if r > p, 



(2.3) 



for some p > 0; p is the radius of the black-hole and is referred to as the "event horizon". 

We now consider, in a coordinate invariant way, Maxwell's equations for an electro- 
magnetic field. First we letf] 

A = Ai dx i 

be a 1-form, called the electromagnetic potential (connection), and the 2- form F = dA 
is the associated electromagnetic field; in coordinates 

F = Fy dx A dx , Fij = -^— - -^— . 



Maxwell's energy momentum tensor is 



Tl3 4vr 



g m Fii>Fj m — -Fi m F m gij 



where as usual, we always use the metric to raise the indices; thus F £m = g h g m ^Fij. 
Maxwell's equations are the pair of equations 

dF = (automatic) 

d*F = , 

where "*" denotes the Hodge star operator, mapping 2-forms into themselves, defined by 



1 
2 



(*F)ki = - ^/\g\e ijM F ij 



where g = det(gij), and e^M is the completely anti-symmetric symbol defined as e^ke = 
sgn(i,j,k,£). It is important to notice that the *-operator depends on the metric (5^). 
The question we now ask is the following: Find the solution of the Einstein/Maxwell 
equations outside of a charged ball in IR 3 . That is, we seek the spherically symmet- 
ric static, exterior (gravitational and electromagnetic) field of a charged distribution of 
matter. This problem was solved in 1918-1920 by Reissner and Nordstrom, and the 
gravitational field is described by the metric 

^fl-^ + ^W-fl-^ + ^r^-rW, 
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Throughout this paper we use summation convention. 



A(r). 




A(r)> 
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Figure 2: The Reissner-Nordstrom Black Hole 

where m = GM/c 2 , q = 2ttGQ 2 /c 4 ; M and Q denote the mass and charge respectively, 
of the ball. Notice that if Q is sufficiently large, the metric coefficient A{r) = 1 — — + -% 
is never zero; otherwise there are two possibilities (see Figure |2|). Case (a) is called a 
non-extreme Reissner-Nordstrom solution and Case (b) is called an extreme Reissner- 
Nordstrom solution; these black-hole solutions will come up again in the next sections. 

We now turn to the Dirac equation. The Dirac equation brings in Quantum Mechanics 
and particles into our work, and it also describes the intrinsic "spin" of the particle 
(fermion). The Dirac equation can be written as 



(G - m) * = 



(2.4) 



where G is the Dirac operator, and \l/ is the wave function of a fermion (proton, anti- 
proton, electron, positron, neutrino, etc.) having (rest) mass m. \l/ is a complex 4-vector 
and is called a spinor. The Dirac operator G takes the form 



G = iGHx) 



d 
dxi 



B(x), 



where the 4x4 matrices G J are called Dirac matrices, and B is a 4 x 4 matrix. The Dirac 
matrices and the Lorentzian metric gi are related by the anti-commutation relations 

gjkj = I| G j 5 G fe| = _(GiG k + G k G j ). (2.5) 

Now let Ti be a space-like hypersurface with a future directed normal vector field v, and 
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define an inner product on solutions of the Dirac equation by 

($, #) = f $G*# i/j d/x , 



m 

where 

is called the adjoint spinor, * denotes complex conjugation, and \x is the invariant measure 
on Ti induced by the metric. This inner product is positive definite and independent of 
Ti because of current conservation, 

V*<K7* = 0. 

From a direct generalization of ^7°v|> in Mf (Minkowski space), where 7° — I n > 

^G^Vj is interpreted as a probability density and for non-black-hole solutions, we nor- 
malize solutions of the Dirac equation by the requirement 

For a black-hole solution with event horizon p, we demand that (cf. [6,7]), 

/ ^G j ^Uj dfi < oo, for all r > p. (2.6) 

J {t=const.,r>ro} 

Now a result in [3] allows us to choose the Dirac matrices G^ to be any 4x4 matrices 
which are Hermitian with respect to the inner product 

iR 4 v 

and satisfy (2.5). Here B(x) is defined by 

B(x) = G j (x)E j (x)+G j (x)A j (x), 
where A^dx 1 is the em potential, 

Ej = \ P(djp) - ^ Tr(G m V 3 G n ) G m G n + % - Tr{ P G 3 V m G m ) P 
l lo o 



is the spin connection, and 



i 



P= : -^\g\e lM G l &G k G". 



In this framework, the Einstein-Dirac-Maxwell (EDM) equations for one particle are 

1 



ij 9 9w ' 


= aT^ 


(G-m)V -- 


= 


V k F> k -- 


= 4rre^G j ^ 



where Tij is the Dirac energy-momentum tensor, and where e denotes the charge of 
the Dirac particles. These are 18 PDE's for the 18 unknowns ^(10), Aj(4), and x I / (4); 
where \P is complex. In this generality, the equations are hopelessly complicated, and are 
impossible to analyze. We thus specialize by imposing certain symmetry conditions; see 
[5, 6] for complete details. 

For j = |, |, .... we consider a static, spherically symmetric system of (2j + 1) Dirac 
particles, each having angular momentum j. (In the language of atomic physics, we 
consider the completely filled shell of states with angular momentum j. Classically, 
one can think of this multi-particle system as several Dirac particles rotating around a 
common center such that their total angular momentum is zero.) Since the system of 
fermions is spherically symmetric, we can obtain a consistent set of equations by assuming 
that both the gravitational and em fields are spherically symmetric. This allows us to 
separate out the angular dependence and reduce the problem to a system of nonlinear 
ODE's. Thus taking the metric in the form (2.2) and the em potential in the Coulomb 
gauge, A = 4>(r) dt, we can take for the Dirac matrices the following: 

G l = T 7 ° 

G r = a/A(7 1 cos i? + 7 2 sin $ cos Lp + 7 s sin "d sin (p) 

I 
G® = -( — 7 1 sin-# + 7 2 cos $ cos 9? + 7 s sin ■# sin 93) 
r 

G^ = -( — 7 2 sin ip + 7 3 cos<^) 

rsinv 
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where 



„_l I \ „ I a' 



T"= _j . V = I, „ , > = 1,2,3, 



• 



and a\a\a are the Pauli matrices (the standard basis for su(2)). 

By using a suitable ansatz and taking the wave functions in the form \1/ = e tuJt f(r 
we can reduce the complex 4-spinors to real two spinors 

$ = (a,/3), a,j3 real. 

The Dirac equations become 

A a' = -^Jl a -[(uo- e0)T + m}(3 (2.7) 

A 13' = \(u- e <i>-m)T-m]a--i-^P, (2.8) 

2r 

and the normalization condition (2.6) is the requirement that 

roc ,/ 7 f ; 

/ (a 2 + (3 2 ) y-=-dr < oo, (2.9) 

for every r > p where r = p defines the event horizon. 

The full EDM equations consist of (2.7), (2.8), together with the Einstein equations 

rA! = 1 - A- 2(2 j + 1) (lu - e0) T 2 (a 2 + (3 2 ) - r 2 AT 2 : (0') 2 (2.10) 

T" f2i-l-T) 2 

IrA— = I- A- 2{2j + l)(u- ecj>) T 2 {a 2 + {3 2 ) + 2 K J - ' T a/3 (2.11) 

+ 2(2j + l)mT(a 2 + p 2 )+r 2 AT 2 ((P') 2 , 
and Maxwell's equation 

r 2 A<f>" = -(2j + 1) e (a 2 + /3 2 ) - UrA + r 2 A ^ + y ^ 0'. (2.12) 

Notice that the EDM equations are invariant under the gauge transformations 

— » + k, uj — > c<j + e/c, fc G H. 
11 



Finally, in addition to (2.9), it is required that the following conditions hold: 

limr(l- A(r)) < oo , (2.13) 

r— >oo 

(finite total mass), 



r— »oo 



Jim T(r) = 1 (2.14) 

(asymptotic flatness), 



lim 0(r) = (2.15) 

(em potential vanishes at infinity). 

We make the following 3 assumptions on the regularity properties of the event horizon 
r = p: 



(I) The volume element v/det \g\ is smooth and non-zero on the horizon: T 2 A 1 and 

T 2 Ae C°°([p, oo)). 

(II) By definition, the strength of the em field is the scalar FijF lj = —2(<J)') 2 AT 2 . We 
assume this to be bounded near r = p, so that in view of (I), we have 

\<p\r)\ < ci, p<r < p + e ± . 

(Ill) A(r) satisfies a power law near r = p, 

A(r)=c(r-p) s + 0(r-p) s+ \ s > 0. 

We remark that if (I) or (II) were violated, then an observer freely falling into the black- 
hole would encounter strong forces when crossing the event horizon; such as an effect is 
seen to be false when on passes to Kruskal coordinates. Hypothesis (III) is a technical 
condition which includes all known physically relevant horizons; s = 1 for Schwarzschild 
and Reissner-Nordstrom black-holes, and s = 2 for the extreme Reissner-Nordstrom 
black-hole. This condition could probably be weakened. 

Here is our main result. 
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Theorem . Every black-hole solution of the EDM equations satisfying Assumptions (I) - 
(III) is either 

a) a non-extreme Reissner- Nordstrom solutions with a(r) = = (5{r), for all r > p, 

or 

b) s = 2 and we find numerically that either the normalization condition (2.9) fails, 
or the solution is not regular everywhere outside the event horizon. 

Thus the only black-hole solutions of the EDM equations are Reissner-Nordstrom 
solutions: the spinors vanish identically, and so are not normalizable. Thus these "quan- 
tum effects" dis-allow (stationary) black-holes. The result indicates too that the Dirac 
particles must either disappear inside the event horizon, or tend to infinity; namely there 
is zero probability for the particles to remain in a finite region outside of the black hole. 

We now give some ideas of the proof; the reader should consult [6] for the full details. 
We assume that the EDM equations admit a solution outside of the even horizon r = p, 
with (a(r),/3(r)) ^ (0,0), and we shall obtain a contradiction. We begin with the case 
s < 2. 

Lemma 1 . If s < 2, then there exist constants c > and e > such that 

c < a 2 {r)+(3 2 {r) <-, p<r<p + e. (2.16) 

c 

Proof . From the Dirac equations (2.7), (2.8) we have 



(^ ~m 



a 



A(a 2 + 2 )' = 2(q,/3) " < l/4m'+ (2jH : 1)2 [> 2 +. 82 ) . 

By uniqueness, (a 2 + (3 2 )(r) > on p < r < p + e for any e > 0. Dividing by VA(a 2 + (3 2 ) 
integrating from r > p to p + e, and noting that A~ x l 2 is integrable near p gives 

, (a 2 + [3 2 )(p + e) "^ / "I-" ' !^ 1/2 

log ■ 



( a 2 + /32)( r ) 



< r p + ^±D! y A -nt)dt 
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< r(^ + ( M±ny /2 A -u Ht)dt<X: 



which implies the result. 

Proposition 2. s < 2 cannot hold. 

Proof . From Einstein's equations (2.10), (2.11), we can write 
r(AT 2 )' = -4(2j + l)(u - e<f))(a 2 + (3 2 )T A 

2 (2j + 1) a(3 + 2(2 j + l)m(a 2 - /3 2 ) 



rpj 



By (I), the lhs is bounded, and T — ► oo as r \ p. Thus, in view of (2.16), we must have 

lim (u - e<f)(r)) = 0. (2.17) 

We next write Maxwell's equation (2.12) in the form 

*• = -i 2 ^ |Q ' 2 + " 2 i - (sb [r2rVlr *") • (2 - 18) 



Since the term ( ) is bounded near r — p, and s > 1 implies A -1 is not integrable, we 
would conclude if s > 1 that 0' is unbounded, thereby violating (II). It follows that we 
must have s < 1. 

Now if we integrate (2.18), we obtain, for r near p, 

0'(r) = Cl (r-p)- s+1 + c 2 + 0((r-p)- s+2 ). 
Integrating again, and using (2.17) gives 

0(r) = Cl (r-p)- s+2 + c 2 (r-p) + - + 0((r-p)' s+3 ) , 

e 

and thus (a; — e<p) = 0(r — p). Now from (2.10) 

rA > = {l-A)-2{2j + l){uj-e<P)[T 2 {a 2 + (3 2 )}-[r 2 {AT 2 ){<P') 2 } 
= O(l) - 0((r - p) 1 -*) - O(l) 
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so that the rhs is bounded, but s < 1 implies that the lhs blows up. This contradiction 
shows that s < 2 cannot hold. 

Proposition 3. s > 2 cannot hold. 

Proof . In [6], we prove the following estimates 

lim(r- p)-i(a 2 (r) + l3 2 (r)) = 0, (2.19) 



and 



lim0'(r) = — lim A _1 (r) T (0 > 0. (2.20) 

r\p p 2 r\p 

Then (2.19) shows that near r = p, 

{ui — e<f>)(r) = c + d (r — p) + o(r — p), d ^ 0, 

and 

[(u - e4)T\' = (r - p)~' [d + c(r - p)- 1 ], 

so that as r \ p, 

[(u — e<f))T] — >■ oo monotonically. 

Now the Dirac equations (2.7), (2.8) can be written as 

a \ ' ( ZLt± -[(u- e0)T + m]\ fa 

A " 

" pj \ [(u - e<f>)T - m] =&±± J \0 

a —b\/a 

b -a J \ 

for r near p, where b /* oo. The eigenvalues of this last matrix satisfy A 2 ~ a 2 — b 2 — *■ — oo 
as r \ p. This suggests that the vector (a, ft) spins around the origin faster and faster 
as r \ p, and that (a, P)(r) -/-> (0, 0). In fact, we prove that 



lim rXp (« 2 (r)+/3 2 (r)j>0, 
but this contradicts (2.19). Thus Proposition 3 holds. 
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The remaining case is s = 2. In this case the metric and electric field behave near 
the horizon like the extreme Reissner-Nordstrom solution. Physically, this means that 
the electric charge of the black-hole is so large that the electric repulsion balances the 
gravitational attraction and prevents the Dirac particles from entering the black hole. 
Obviously, this is not the physical situation which one can expect, say, from the gravi- 
tational collapse of a star. On the other hand, extreme Reissner-Nordstrom black-holes 
have zero temperature [10, 17], and can be considered as end-states of black-holes emit- 
ting Hawking radiation. It is thus interesting to see whether the case s = 2 can admit 
normalizable solutions of the EDM equations. Considering the Dirac equation in an 
extreme Reissner-Nordstrom black-hole background field, it was proved in [6] that the 
normalization condition (2.6) is violated. What we want to know here is whether the 
influence of the spinors on the gravitational and em fields can make the normalization 
condition finite. Clearly this is a very difficult mathematical problem. We have made 
numerical investigations, and from these we conclude that the solutions of the EDM equa- 
tions either develop singularities at some finite r, or the normalization condition fails. 
Thus in the case s — 2, our numerical investigations show that there are no normalizable 
solutions of the EDM equations. 

Our results show that taking quantization and spin into account implies a breakdown 
of the classical situation; namely, there cannot exist normalizable black-hole solutions 
of the EDM equations. Applied to the gravitational collapse of a "cloud" of spin-1/2 - 
particles to a black hole, our result indicates that the Dirac particles must either disappear 
into the black-hole, or escape to infinity. 

We remark that these results have been extended to the axisymmetric case including 
the Kerr-Newman rotating black hole [9]. In addition, we have recently proved that 
there are no normalizable solutions of the Einstein-Dirac-S , ?7(2) Yang/Mills equations; 
this result holds under a more general condition than (III), and applies for every s > 0. 

We close with two facts. First, our results are basically a consequence of the Heisen- 
berg Uncertainty Principle, together with the specific form of the Dirac current. Second, 
it is essential for our results that the particles have spin; spin is an important effect which 
must be considered in studying gravitational collapse. 
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